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Lattice Kirigami, ultra-light metamaterials, poly-disperse aggregates, 
ceramic nano-lattices, and two-dimensional (2-D) atomic materials 
share an inherent structural discreteness, and their material prop¬ 
erties evolve with their shape. To exemplify the intimate relation 
among material properties and the local geometry, we explore the 
properties of phosphorene -a new 2-D atomic material- in a conical 
structure, and document a decrease of the semiconducting gap that 
is directly linked to its non-planar shape. This geometrical effect oc¬ 
curs regardless of phosphorene allotrope considered, and it provides 
a unique optical vehicle to single out local structural defects on this 
2-D material. We also classify other 2-D atomic materials in terms 
of their crystalline unit cells, and propose means to obtain the lo¬ 
cal geometry directly from their diverse two-dimensional structures 
while bypassing common descriptions of shape that are based from 
a parametric continuum. 

Phosphorene | Discrete differential geometry | Two-dimensional materials 
Abbreviations: DDG, discrete differential geometry; 2-D, two-dimensional 

T wo-dimensional (2-D) materials [1-20] are discrete sur¬ 
faces that are embedded on a three-dimensional space. 
Graphene [1, 2] develops an effective Dirac-like dispersion on 
the sublattice degree of freedom and other 2-D atomic materi¬ 
als exhibit remarkable plasmonic, polariton, and spin behav¬ 
iors too [18-20]. 

The properties of 2-D materials are influenced by their lo¬ 
cal geometry [12-17,21-29], making a discussion of the shape of 
two-dimensional lattices a timely and fundamental endeavor 
[24, 30, 31, 32]. A dedicated discussion of the shape of 2-D 
materials is given here within the context of nets. Nets are 
discrete surfaces made from vertices and edges, with vertices 
given by particle/atomic positions [31-34]. The discrete ge¬ 
ometry that originates from these material nets is richer than 
its smooth counterpart because the net preserves the struc¬ 
tural information of the 2-D lattice completely, yielding exact 
descriptions of shape that remain accurate as the lattice is 
subject to arbitrarily large structural deformations [15, 16], as 
the particle/atomic lattice becomes the net itself. 

A discrete geometry helps addressing how strain influences 
Chemistry [35], how energy landscapes [36, 37, 38] correlate 
to non-planar shapes, and it provided the basis for a lat¬ 
tice gauge theory for effective Dirac fermions on deformed 
graphene [39, 40]. In continuing with this program, the optical 
and electronic properties of phosphorene cones will be linked 
to their geometry in the present work. At variance with 2- 
D crystalline soft materials that acquire topological defects 
while conforming to non-planar shapes [12-15], the materials 
considered here have strong chemical bonds that inhibit plas¬ 
tic deformations for strain larger than ten percent [41]. 

The study is structured as follows: We build conical struc¬ 
tures of black and blue phosphorene, determine their local 
shape and link this shape to the magnitude of their semicon¬ 



Fig. 1 . Unit cells for (A) black and (B) blue phosphorene monolayers that are 
formed by two sublayers (Si and S 2 ) separated by a distance To. Sublayers in black 
phosphorene are irregular -distances among nearest atoms belonging to the same sub¬ 
layer are not identical- as emphasized by the dashed lines joining atoms belonging to 
sublayer Si on the rightmost structure in (A). 

ducting gap. It is clear that a discrete geometry applies for 
arbitrary 2-D materials. 


Results and discussion 

Phosphorene [8, 9] has many allotropes that are either semi¬ 
conducting or metallic depending on their two-dimensional 
atomistic structure [42, 43, 44]. The most studied phase, 
black phosphorene (Fig. 1A), has a semiconducting gap that 
is tunable with the number of layers, and by in-plane strain 
[26, 27, 45, 46]. Theoretical studies of defects on planar phos¬ 
phorene indicate that dislocation lines do not induce localized 
electronic states [47], and algorithms to tile arbitrary planar 2- 
D phosphorene patterns have been proposed as well [44]. Unit 


Significance 

Phosphorene is a new 2-D atomic material, and we document 
a drastic reduction of its electronic gap when under a conical 
shape. Furthermore, geometry determines the properties of 2- 
D materials, and we introduce discrete differential geometry to 
study them. This geometry arises from particle/atomic posi¬ 
tions; it is not based on a parametric continuum, and it applies 
across broad disciplinary lines. 
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Fig. 2. (A) Black and (B) blue phosphorene cones are built by removing the angular segments in white from the planar structures (as illustrated at the upper-right corner 

of the conical structures), joining the cut structures along the red (disclination) lines, and a subsequent structural optimization. The discrete local geometry of these conical 
structures at sublayers Si and S 2 is given by four local invariants ( Tr(g ), Det(g), H and K) that are obtained at each atomic position. These invariants are contrasted 
with the geometry of planar structures ( Tr(g ) = 1, Det(g) = 1, H = 0, and K = 0) that is depicted at the uppermost row for sublayer S±. The change in relative 
height t/to is shown at each atomic position as well. Features such as: (i) strain induced by the dislocation line in planar black phosphorene; (ii) a sensible compressive strain 
near the apex; (iii) curvature lacking a radial symmetry; and (iv) the lack of significant changes in r on conical black phosphorene, are clearly visible. The disclination line 
on the blue phosphorene cone is reflected on the metric invariants and on H\ This disclination has a semi-cylindrical shape, that yields an overall radially-symmetric Gaussian 
curvature K. In addition, a decrease of r, and an in-plane compression (elongation) at sublayer Si (S 2 ) near the apex can also be seen on the blue phosphorene cone. 


cells of planar monolayers of black and blue phosphorene are 
displayed on Fig. 1. 

Four invariants from the metric ( g ) and curvature ( k ) ten¬ 
sors determine the local geometry of a 2-D manifold [24, 30]: 

Tr(ff), Det(g), = Tr(fc)/2Tr(<?), K = Det(fc)/2Det(<?), [1] 

where Tr (Det) stands for trace (determinant), H is the mean 
curvature, and K is the Gaussian curvature. 

An infinite crystal can be built from these unit cells, and 
the geometry of such ideal planar structure can be described 
by Tr(g) = 1, Det(g) = 1 (i.e., no strain), H — 0, and K — 0 
(i.e., no curvature) at both sublayers Si and S 2 . In addition, 
their thickness r is equal to To before a structural distortion 
sets in: These five numbers quantify the local reference, flat 


geometry. (In principle, Tr(g) and Det(g) will be functions 
of the lattice constant, but here they are normalized with re¬ 
spect to their values in the reference crystalline structures seen 
on Fig. 1, in order to enable direct comparisons of the met¬ 
ric among black and blue phosphorene monolayers. To = 2.27 
(1.26) A for black (blue) phosphorene in these planar reference 
structures.) 

Now, structural defects can induce non-zero curvature and 
strain [35], and may also be the culprits for the chemical degra¬ 
dation of layered phosphorene. In order to study the con¬ 
sequences of shape on the optical and electronic properties 
of phosphorene monolayers, we create finite-size conical black 
and blue phosphorene monolayers, characterize the atomistic 
geometry, and investigate the influence of shape on their semi¬ 
conducting gap. (We have proposed that hexagonal boron 
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nitride could slow the degradation process while allowing for 
a local characterization of phosphorene allotropes [48], and a 
study of chemical reactivity of non-planar phosphorene will be 
presented elsewhere.) 

Phosphorene cones are built from finite disk-like planar 
structures that have hydrogen-passivated edges (c.f., Fig. 2). 
The black phosphorene cone seen in Fig. 2A is created as fol¬ 
lows: We remove an angular segment -that subtends a 4> = 46° 
angle- from a planar structure that has a dislocation line [47]. 
The two “ridges” that are highlighted by the red segments 
on the planar structure in Fig. 2A are joined afterwards to 
create a disclination line. Atoms are placed in positions dic¬ 
tated by an initial (analytical) conical structure, and there 
is a full structural optimization via molecular dynamics at 
the ab-initio level to relieve structural forces throughout (see 
Methods). 

Blue phosphorene has a (buckled) honeycomb structure 
reminiscent of graphene, so the conical structure seen to the 
left of Fig. 2B was generated by removing an angular segment 
subtending a <f> = 60° angle on planar blue phosphorene, and 
following prescriptions similar to those employed in creating 
the black phosphorene cone afterward. 

The subplots displayed to the right on Figs. 2A and 2B show 
the local discrete geometry at individual atoms (c.f., Eqn. [1]; 
see details in the Methods section). For each allotrope, the 
data is arranged into three rows that indicate the geometry of 
the planar structure at sublayer Si , and the local geometry of 
the cones at sublayers Si and S 2 • An additional plot shows 
the value of t/to that tell us of local vertical compression. 

There is strain induced by the dislocation line on the pla¬ 
nar black phosphorene structure, as indicated by the color 
variation on the Tr(g) and Det(g) plots, that implies having 
atoms at closer/longer distances than in the reference struc¬ 
ture, Fig. 1. A slight curvature on the black phosphorene 
planar structure, induced by the dislocation line, is also vis¬ 
ible on the H plot in Fig. 2A. The planar blue phosphorene 
sample does not have any dislocation line, and for that reason 
the metric shows zero strain ( Tr(g ) = 1 and Det(g) = 1) and 
zero curvature (H = 0 and K = 0) on that reference structure 
(Fig. 2B). 

The black phosphorene conical structure seen on Fig. 2A 
has the following features: A compression near its apex, as dis¬ 
played by the white color on the metric invariants; this com¬ 
pression is not radial-symmetric. An asymmetric elongation 
is visible towards the edges. In addition, there is a radially- 
asymmetric non-zero curvature; the observed asymmetry re¬ 
flects the presence of the dislocation/disclination axis, which 
provides the structure with an enhanced structural rigidity. 
This rigidity is confirmed by the ratio r/ro close to unity, 
which indicates almost no vertical compression/elongation of 
this conical structure (Fig. 2A). 

The blue phosphorene cone has a more apparent radial sym¬ 
metry, except for the disclination line that is created by the 
conical structure, as reflected on the metric invariants and on 
H. There is compression (elongation) at the apex, and elon¬ 
gation (compression) along the disclination line in sublayer Si 
(S' 2 ). The disclination line has a semi-cylindrical shape and 
hence a zero radius of curvature along the disclination line, 
resulting in a zero Gaussian curvature along such line; the 
Gaussian curvature looks radially-symmetric overall. 

Figure 2 indicates that the blue phosphorene cone acquires 
the largest curvatures of these two cones. This is so because 
the angular segment removed from the planar blue phospho¬ 
rene sample has a comparatively larger value of 0. One notes 



Fig. 3. The semiconducting gap of planar black and blue phosphorene monolayer 
discs increases as the number of atoms decreases. Even though the number of atoms 
must be reduced to create a conical shape from a planar structure, the gap decreases 
on conical structures with such decrease of the number of atoms, as indicated by open 
ticks and the tilted arrows: This reduction of the gap on a conical structure thus has 
a purely geometrical origin. (Dash and dash-dot lines indicate these gaps when the 
number of atoms is infinite.) 

the rather smooth curvature at the apex on the blue phospho¬ 
rene conical structure after the structural optimization. 

The change of r with respect to To is created by an out-of¬ 
plane strain or by shear. The blue phosphorene cone shows 
out-of-plane compression near the apex. The distance among 
planes in black phosphorene is closer to its value in an ideal 
planar structure throughout, showing scatter around the dis¬ 
location line. 

The main point from Figure 2 is the strain induced by curva¬ 
ture. The strain pattern observed on that Figure is far more 
complex than those reported before for planar phosphorene 
[26, 27, 45, 46], and the discrete geometry captures the strain 
pattern with the precision given by actual atomic positions. 
We will describe the tools that lead to this geometry based 
from atoms later on (Methods), but we describe the effect of 
shape on the material properties of these cones first. 

Black and blue phosphorene monolayers are both semicon¬ 
ducting 2-D materials with a direct bandgap, and we inves¬ 
tigate how this semiconducting gap evolves with their shape. 
The semiconducting gap E g — is equal to 1.1 (2.2) eV for the 
finite-size planar black (blue) phosphorene monolayers on dis¬ 
play in Fig. 2, before the removal of the angular segments. 
These electronic gaps are highlighted at the rightmost end of 
Fig. 3. All gaps were computed after a structural optimiza¬ 
tion, in order to avoid unbalanced forces on these samples that 
would bias their magnitude (Methods). 

We determine size-effects on the semiconducting gap on pla¬ 
nar structures first: The magnitude of the gap E g increases as 
the number of atoms decreases on strain-free disk-like planar 
structures; c.f., Fig. 3. For an infinite number of atoms —i.e., 
for a fully periodic planar 2-D crystal- the gap converges to 
the values indicated by the dash and dash-dot lines on Fig. 3; 
namely, 0.8 eV and 2.0 eV for black and blue phosphorene, 
respectively. These magnitudes were obtained from standard 
density-functional theory [8, 42]. Although we acknowledge 
that other methods describe the dielectric properties more ac¬ 
curately on smaller samples [27], our focus is on the trend the 
gap follows, and given that the trend is geometrical in nature, 
it will stand correct despite of the particular method employed 
in computing the semiconducting gap. 

The conical structures have less atoms than their parent pla¬ 
nar structures once the angular segments seen in Fig. 2A and 
2B are removed. We learned on previous paragraph that the 
electronic gap increases as the number of atoms is decreased 
on planar structures. Following that argument, one may ex- 
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Fig. 4. A remarkable reduction of the semiconducting gap occurs due to the strain induced by the conical shape, and regardless of phosphorene allotrope. The local 
distribution of A Ei correlates with the extent of the electronic wavefunctions shown below and above the Fermi level for (A) black and (B) blue phosphorene. 


pect the semiconducting gap for the conical structures to be 
larger than the one observed on their parent planar structures. 

But instead of increasing with the removal of atoms, the gap 
decreases dramatically on the conical structures, going down 
to 0.84 (1.70) eV for the black (blue) phosphorene cone, as 
indicated by the open red symbols, and the tilted arrows in 
Fig. 3: Shape influences phosphorene’s material properties. A 
practical consequence of this result is that optical probes could 
provide a measure of the local shape of phosphorene. The re¬ 
lation among the semiconducting gap and a non-planar shape 
has been observed in transition metal dichalcogenide monolay¬ 
ers recently [25], and we establish here that structural defects 
on semiconducting 2-D materials create a similar effect. Let 
us next address the reason for such reduction of the semicon¬ 
ducting gap. 

The semiconducting gap is a global material property. On 
finite samples it is possible to project the electronic density 
onto individual atoms to learn about the spatial localization 
of the electronic states below and above the Fermi level, whose 
energy difference is equal to the electronic gap. Those states 
are known as the highest occupied molecular orbital (HOMO) 
and the lowest unoccupied molecular orbital (LUMO), and 
they were spin-degenerate in all the samples we studied. Addi¬ 
tionally, the n—th electronic orbital below (above) the HOMO 
(LUMO) is labeled HOMO-n (LUMO+n). 

The HOMO and LUMO states cover the entire planar black 
and blue phosphorene structures (c.f., Fig. 4) and are thus 
delocalized. On the other hand, the conical samples have or¬ 
bitals below and above the Fermi energy that display a certain 
amount of localization. 

Such effect is most visible for the HOMO and HOMO—1 
states for the blue phosphorene cone in Fig. 4B. Thus, unlike 
a dislocation line [47] the pentagon defect that is responsible 
for the curvature of the blue conical structure is localizing the 
HOMO state within a ~ 10 A radius. Given that the pentagon 
defect originates curvature itself, we conclude that curvature 
leads to a reduction of the electronic gap on this system. 

Given the localization observed on some of these orbitals, 
we can further ask: At a given atomic position, what is the 
first orbital with a non-zero density at such location below 


and above the Fermi energy? This question can be rephrased 
in terms of the difference in energy among the first orbitals 
below and above the Fermi level that have a non-zero proba¬ 
bility density at a given atomic location: On the third column 
on Figs. 4 A and 4B we display such energy difference A Ei 
at every phosphorus atom i. There exists a clear correlation 
among A Ei and the localization pattern of the orbitals below 
and above the Fermi energy, as it should be the case. The 
reduction on the semiconducting gap on the cones is empha¬ 
sized by normalizing A Ei in terms of E g , which yields an 
twenty-percent reduction of the semiconducting gap for both 
allotropes: A non-planar geometry that is created by struc¬ 
tural defects on any phosphorene allotrope will lead to a size¬ 
able reduction of the semiconducting gap. 

To end this work, we must state that the discrete geometry 
used to tell the shape of the phosphorene conical structures 
applies to other 2-D atomic materials with varied unit cells; 
some of which are listed below: 

1. Regular honeycomb lattices (graphene [2, 39, 40, 49] and 
hexagonal boron nitride (hBN) [4]. 

2. Low-buckled honeycomb hexagonal lattices (silicene and 
germanene [5]. We established that freestanding stanene 

[50] is not the structural ground state. 

3. “High-buckled” hexagonal close-packed bilayers of bismuth 

[51] , tin and lead [52]. 

4. Thin trigonal-prismatic transition metal dichalcogenides 
[53, 54]. 

5. Materials with buckled square unit cells -quad-graphs [31]— 
such as A1P [7]. 

Low-buckled structures (silicene, germanene, blue phospho¬ 
rus), hexagonal close-packed bilayers (bismuth, tin and lead) 
and black phosphorus have two parallel sublayers Si and S 2 . 
Transition-metal dichalcogenides and A1P have three parallel 
sublayers Si — S3. Structures 1—4 are equilateral triangu¬ 
lar lattices with a basis, for which individual planes represent 
regular equilateral triangular nets; and structure 5 realizes a 
regular quad-graph. The discrete geometry clearly stands for 
other 2-D materials that do not form strong directional bonds 
[12, 14] as well. 
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Concluding remarks 

Shape is a fundamental handle to tune the properties of 2-D 
materials, and the discrete geometry provides the most accu¬ 
rate description of two-dimensional material nets. This ge¬ 
ometry was showcased on non-planar phosphorene allotropes 
for which the electronic gap decreases by twenty-percent with 
respect to its value on planar structures. The discrete ge¬ 
ometry can be thus be used to correlate large structural de¬ 
formations to intended functionalities on 2-D materials with 
arbitrary shapes. 

Methods 

Creation of conical structures. Consider structures having 
atomic positions at planar discs; these positions can be pa¬ 
rameterized by r*i, 0 i , and Zi — {0,to}. Calling the angu¬ 
lar segments being removed from these planar structures, the 
range of the angular variable 0i is [0, 27r — 0], and the conical 
structures are initially built by the following transformation: 

r 'i = n\J 1 - 0'i = 2 and z[ = Zi-n-^. These con- 

ical structures containing about five hundred atoms undergo a 
structural optimization via ab-initio (Car-Parinello) molecular 
dynamics [55] with the SIESTA code [56, 57], until forces are 
smaller than 0.04 eV/A. 

Calculation of the electronic gap of phosphorene samples. We 

obtain A Ei as follows: Let Ep be the charge neutrality level 
or Fermi energy, and pi(E) the density of electronic states 
projected onto atom i. We call Ef (Ef) the first energy level 
observed on pi(E) lying above (below) Ef at atom i, and re¬ 
port A Ei/Eg = (Ef — Ef)/E g on Fig. 4. We note that all 
structures we worked with had a final net zero spin polariza¬ 
tion on spin-polarized calculations. 

A discrete geometry based on triangulations at atomic posi¬ 
tions. There is an extrinsic and continuum (Euclidean) geom¬ 
etry in which material objects exist. These material objects 
are made out of atoms that take specific locations to generate 
their own intrinsic shape. 



Fig. 5. Discrete tensors based on triangulations are expressed in terms of averaged 
normals hj, edges Gj , the normal of the triangle is, and dual edges e * = Gj X is 
(j = 1,2,3). These tensors can be employed in quad-nets, irregular lattices, and 
structures containing defects as highlighted by the pentagon surrounding atom 1. 

The intrinsic shape of 2-D atomic materials can be idealized 
as a continuum: The basic assumption of continuum mechan¬ 


ics is that a continuum shape is justified at length scales l much 
larger than interatomic distances that are characterized by a 
lattice parameter ao ( l >> ao). The common understanding 
of shape arises from within this continuum perspective that is 
based on the differential geometry of two-dimensional mani¬ 
folds. This continuum approximation is valid down to inter¬ 
atomic scales (Z < ao) for slowly-varying deformations, but not 
when curvature concentrates at atomic-scale bonds or pleats, 
such as in the examples provided in Ref. [11] and [14]. And so, 
while these sharp structures may be disregarded or approxi¬ 
mated within the context of a continuum surface, we hold the 
opinion that an intrinsic geometry that is exact at the atomic 
scale must form part of the theoretical tool-set to deal with 
2-D materials. The discrete geometry to be described in the 
next paragraphs bypasses descriptions of shape that are based 
on an effective continuum, and it brings an understanding of 
the shape of material nets at a fundamental level. 

Consider three directed edges ei, e 2 and e 3 , such that 
ei +e 2 +e 3 = 0 , and define Q{ = e/ • e/ (Z = 1,2,3), represent¬ 
ing the square of the smallest finite distance among atoms on 
the 2-D lattice. This is the discrete analog of the infinitesimal 
length element ds 2 (c.f., Fig. 5). 

Now consider the change in orientation among normals hj 
and hfc, and project such variation onto their common edge 
ej: This is is, define Q\ t = (fife — hj) • ei (see Fig. 5; j, /c, l are 
permutations of integers 1, 2, and 3). The reader may recall 
that the curvature tensor is defined as h • , but the discrete 

d£/3 ’ 

tensor carries a constant edge, and it rests on changes of the 
local normals instead. In previous Equation, rp is the average 
over individual normals at triangulated area elements within 
the polygon surrounding atom l and highlighted by dashed 
lines. The dual edge is defined by e* = ei x is , with is the 
normal to the triangle formed by atoms 1 , 2 and 3 and At is 
the triangle area (—e 2 xei = 2 Atv)- 

This way, the discrete metric tensor takes the following form 
on the basis given by dual edges [58]: 

Y (Qj-Qk [2] 

0 0,M) 

with Ao the area of the triangulated area element at the ref¬ 
erence (non-deformed and defect-free) plane. 

The discrete curvature tensor has an identical structure: 

= (QV - Tk 1 - QV)e* « e*. [ 3 ] 

T (. 

The parenthesis (j, k,l) indicates a sum of three terms, as fol¬ 
lows: (j = 1, k = 2, l = 3), (2, 3, 1), and (3, 1, 2). Eqns. [2] 
and [ 3 ] become 3x3 matrices with explicit values for Qj , Qj 1 
and e* from atomic positions (Fig. 5). For instance, the dis¬ 
crete curvature tensor has eigenvalues {0, k±, Zc 2 } at each trian¬ 
gulated area element, yielding H = (Aq + fo)/2 and K = kik 2 . 
The geometrical invariants reported at point j are averages 
over their values at individual triangles sharing this vertex. 

Topological defects (exemplified by a pentagon seen by a 
dot-dash line on Fig. 5) break translation symmetry, making 
it impossible to recover a crystalline structure by means of 
elastic deformations. Yet, Eqns. [ 2 ] and [ 3 ] provide a geom¬ 
etry along topological defects seamlessly still, as seen in the 
geometry provided in Fig. 2. 

The computation of the distribution of the thickness r in 
both conical structures was performed by finding for each atom 
in the lower sublayer Si the three nearest neighbors on the up¬ 
per sublayer S 2 . The distance between the centroid these three 
atoms and the atom in question at the Si sublayer amounts 
to the thickness r of the structure at that specific atom. 
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